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SECTION  I 
INTRODUCTION 


In  several  previous  publications  by  the  Numerical  Analysis  Group  at  Swansea  (References 
If  2,  3,  4,  and  5)  the  methodology  of  generating  curved,  ‘  isoparametric’  elements  was  fully 
explained.  By  suitable  transformations,  curvilinear  coordinates  within  such  elements  are 
generated  and,  using  numerical  integration,  the  element  properties  can  be  determined. 
Two-dimensional  problems  of  plane  stress  or  strain,  axisymmetric  problems  (Reference  3), 
and  fully  three-dimensional  problems  (References  4  and  5)  can  be  treated,  using  elements 
of  the  families  shown  in  Figure  1  and  2. 

Evidently  curved  shells  can  be  attempted  simply  by  making  such  two-  or  three-dimensional 
elements  thin.  Indeed  up  to  a  point  this  has  proved  successful  in  practice;  for  instance,  an 
axisymmetric  solid  program  based  on  the  elements  of  Figure  1,  b  and  c  has  been  used  with 
a  length- thickness  ratio  of  50  or  more.  Two  limitations  arose  however.  In  the  first  place 
round-off  errors  occur  with  very  large  length-thickness  ratios.  In  the  second  place,  an 
unnecessarily  large  number  of  degrees  of  freedom  are  carried,  because  simplifying  con¬ 
straints  are  usually  introduced  in  shell  analyses.  For  instance,  even  in  an  axisymmetric 
problem  using  elements  of  Type  1c,  there  are  8  degrees  of  freedom  within  the  shell  thickness, 
whereas  with  thin  shell  assumptions  only  three  would  strictly  be  necessary.  Further,  for 
membrane  situations  the  number  of  degrees  of  freedom  may  be  reduced  even  more. 

This  paper  deals  with  degenerate  versions  of  these  general  elements  for  shell  problems. 
In  the  early  sections  axisymmetric  shell  situations  are  treated  and  later  the  process  is 
extended  to  cover  general  shells. 

In  the  ‘degeneration’  process,  two  assumptions  are  made.  Firstly,  the  original  normals 
to  the  middle  surface  are  taken  as  inextensible  and  straight.  Secondly,  the  elastic  modulus 
in  the  normal  direction  is  taken  as  zero,  to  avoid  normal  stresses. 

The  above  assumptions  are  compatible  with  thin  shell  theory,  but  they  allow  the  normals 
to  rotate  relative  to  the  middle  surface  after  deformation.  Thus  a  measure  of  shear  deforma¬ 
tion  of  the  shell  is  retained.  Introducing  a  slight  ‘  correction  factor’  to  the  shear  deformation, 
it  appears  that  the  answers  are  consistent  with  Reissner’s  (Reference  6)  theory  of  plate  in 
practical  applications. 
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Parent  - Generalized 


Figure  1.  Two  Dimensional,  Quadrilateral,  Family 
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These  new  elements  are  not  strictly  ‘isoparametric’  according  to  the  definitions  of 
Reference  3.  Nevertheless  conformity  between  elements  is  maintained,  and  rigid  body  modes 
and  constant  strain  criteria  are  in  turn  exactly  and  approximately  satisfied. 


SECTION  n 

COORDINATE  SYSTEM  FOR  AN  AXISYMMETRIC  ELEMENT 


Consider  a  system  of  curvilinear  coordinates  £  -  tj  defined  inside  the  element  shown  in 
Figure  3a,  with  typical  node  pairs  i^  and  bottom'  The ‘normal’  direction  rj  as  defined  by 
such  points  is  in  general  only  approximately  normal  to  the  mid-surface.  However,  it  is  in 
this  direction  rj  that  the  linearity  constraint  is  imposed. 

The  coordinates  r  and  z  within  the  element  are  linked  with  the  curvilinear  coordinates 
by  the  relation 


r 

=  In.  ( 1  +  1?) 

r. 

i 

+  1",  ^ 

~ri  ' 

z 

1  2 

Zj 

(i) 

- 

_ 

top 

L  ■  J 

bottom 

Where  N.  is  a  shape  function  in  £  only.  The  element  extends  from  £  and  f]  =  -1  to  1  in  the 
usual  manner  (Reference  3). 


This  can  also  be  written  in  terms  of  the  mid-surface  coordinates  and  the  angle  of  the 
‘  normal’  shown  in  Figure  3b  as  follows: 


r 

=  I  N 

ri 

r*  fi 

+  I  N  7f  -L 

cos  ^>. 

2 

1 

_Zi_ 

1  '  2 

mid 

s  in  <£. 

where  t.  is  the  scalar  length  of  the  normal. 
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It  is  convenient  to  use  linear  functions  N.  for  elements  defined  by  two  nodes,  parabolic 
functions  for  three  nodes  and  so  on  as  shown  in  Figure  4,  so  that  middle  surfaces  of  elements 
become  linear,  parabolic,  cubic.  Typical  shape  functions  are 


Linear 


Line  ar 


Where  £ 


i  ■  T  (l+«o> 

CCi  and  4 •  =  ±  I 


(3) 


Quadratic  For  an  end  node 

Ni  =  \  (  1  +Zo] 

and  for  a  mid  side  node 

n  j  =  d-C2) 

Cubic  For  an  end  node 

Ni  =  (  n-C  >  t-i  +  90 

1  16  o 

and  for  a  node  with  £  =  ±  -i- 

Nj  =  tl  -£*  )  (I  + 9  Co  > 


(4) 

(5) 

(6) 

(7) 
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SECTION  m 

DISPLACEMENT  PATTERN  OF  AN  AXISYMMETRIC 
ELEMENT  UNDER  AXISYMMETRIC  LOADING 


Assuming  that  the  ‘  normals’  are  straight  and  non- extensible,  the  general  displacements 
(u,  v)  are  defined  by  the  displacements  of  the  mid-surface  nodes  in  the  r  -  z  directions  and 
the  rotation  of  the  normal  at  such  nodes  (u.,  v.  and  cr  of  Figure  5).  Within  the  element  the 
displacements  thus  defined  are 


u 

U. 

1 

’  -  Sin 

V 

•In, 

.  vi . 

- 

+  I  Nj  T) 

COS  <£j 

(8) 


Comparison  with  Expression  2  shows  that  the  isoparametric  concept  is  retained  only  in 
the  £  direction. 


Figure  5.  'Global*  Displacements  in  an  Axisymmetric  Shell 
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SECTION  IV 

EVALUATION  OF  ELEMENT  CHARACTERISTICS 
FOR  AXISYMMETRIC  PROBLEMS 


Element  characteristics  such  as  stiffness  matrices  etc.,  can  now  be  determined  in  the 
usual  manner.  For  example,  the  stiffness  matrix  is  defined  as 

[kf  ./  [b]T  [d]  [b]  d(Vol)  (9) 

Vol. 

where  the  {^B j  matrix  relates  the  local  strains  to  the  nodal  displacements 

W  ■  [■]{»}' 

using  the  notation  of  the  finite  element  text  (Reference  7). 


However,  it  is  necessary  to  avoid  the  spurious  side-effects  of  inextensibility  by  assuming 
a  zero  modulus  in  the  true  normal  direction.  The  local  strains  must  therefore  be  defined 
in  appropriate  directions. 


Consider  local  orthogonal  coordinates,  z/  parallel  to  a  surface  V  =  constant  within  the 
shell,  and  r'  created  truly  normal  to  this,  with  the  corresponding  displacement  components 
v'and  u'  as  shown  in  Figure  6,  The  strains  in  these  local  coordinates  are  given  by 

111 

dz' 

u 

r 

d  v'  _du ' 
dr'  dz' 

neglecting  the  strain  normal  to  the  surface  in  accordance  with  our  shell  assumptions,  and 
writing  the  hoop  strain  directly  in  global  coordinates. 

It  is  now  necessary  to  write  these  strains  in  terms  of  the  nodal  displacements  and  the 
local  coordinates  £  ,  .  By  Equations  2  and  8  the  second  term  is  directly  available;  but 
two  transformations  are  needed  to  obtain  the  others.  These  are  easily  derived. 
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Figure  6.  Local  Orthogonal  Coordinates  for  Symmetric  Loading 


Firstly,  the  r,  z  derivatives  of  u,  v  are  derived  from  their  17  derivatives.  This  pro¬ 
cess  is  well  known  and  gives 
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with  the  Jacobian  matrix  J  defined  as 


d  z 
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ae 

ae 

dz 

dr 
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(12) 


(13) 


Secondly,  changing  to  local  orthogonal  coordinates  one  can  write 
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where  [0J  is  the  matrix  of  direction  cosines  of  the  z'  and  r'  axes.  Parametrically  these  can 
be  obtained  at  any  value  of  £  and  77,  noting  that  directions  z'  and  £  coincide,  as 


In  thi s_,  all  the  derivatives  are  obtained  numerically  from  Equation  2. 

The  |^B']  matrix,  which  defines  the  local  strains  in  terms  of  the  nodal  displacements,  can 
now  be  evaluated  at  any  point  within  the  element. 


{«'}  ■  [■'] 


a, 


CX; 


[»■]  {*}! 


(16) 


The  corresponding  elasticity  matrix  jp'j  is  now  easily  written.  In  an  isotropic  case 
for  instance  it  is 


M 


i 

V 


V 

I 


0 

0 

I  -  v 
2k 


(17) 


Where  E  and  v  are  the  elastic  modulus  and  Poisson’ s  ratio  respectively.  The  factor  k  is 
included  to  account  more  accurately  for  shear  strain  energy.  As  the  displacements  vary 
linearly  across  the  thickness  of  the  shell,  the  shear  stresses  are  sensibly  constant.  However, 
they  are  known  to  be  approximately  parabolic,  and  the  factor  k  =  1.2  is  included  to  improve 
the  representation  (see  standard  strength  of  materials  theory  of  shear  deformation  energy). 


The  final  step  is  to  define  the  volume  element  for  the  integrations,  such  as  that  in  the 
stiffness  matrix  of  Equation  9.  This  again  is  standard;  for  an  annular  element  it  becomes 

d(vol)=27rr(detj-jjd£d77)  (|8> 
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The  actual  integration  is  carried  out  numerically  within  the  simple  limits  of  ±1,  using 
several  Gauss  points  as  described  in  References  1  to  3.  Only  two  Gauss  points  are  needed 
in  the  17  (thickness)  direction,  and  3  or  4  in  the  £  direction,  depending  on  whether  the  elements 
are  parabolic  or  cubic. 

Either  Equation  1  or  2  can  be  used  to  define  the  element  geometry  in  the  program.  For 
very  thick  shells  Equation  1  may  be  preferable,  but  for  thinner  shells,  particularly  if  branch¬ 
ing  occurs,  the  thickness  slope  definition  is  better. 


SECTION  V 

EXTENSION  TO  NON-SYMMETRIC  LOADING 


It  is  now  well  known  that  nonsymmetric  loads  acting  on  axisymmetric  bodies  can  be 
treated  separately  for  each  harmonic,  using  Fourier  analysis,  because  their  effects  are 
uncoupled  (References  8,  9,  and  10).  Thus,  if  the  radial,  axial  and  tangential  loads  are  given 
by 

R 
Z 
T 


Rn  cos  n 9 

Zn  cos  n  Q 

Tn  sin  n  9 


(19) 


then  the  resulting  displacements  are  given  by 


u 

u  n 

cos 

n0  ' 

V 

= 

vn 

cos 

n0 

w 

"n 

sin 

n$ 

* 

— 

By  considering  the  amplitudes  of  the  displacement,  the  appropriate  stiffness  matrices  can 
be  found  by  a  volume  integration  which  includes  the  Q  variations. 
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To  extend  the  axisymmetric  analysis  in  this  way,  we  must  add  one  further  degree 
of  freedom  w  to  the  general  displacement  in  Equation  8,  and  two  further  degrees  of  freedom 
at  each  nodal  normal,  (one  displacement  and  one  rotation)  giving  five  altogether. 


Assuming  the  loading  is  as  given  by  Equation  19  and  retaining  the  assumptions  used 
previously  in  deriving  the  displacement  pattern,  we  can  write  in  place  of  Equation  8 


u 

cos 

nd  0 

0  1 

I 

u. 

1 

-sin  <£. 

0" 

a. 

\ 
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0 

cosn£? 
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+ 1  Ni  V  T 

COS  <£j 

0 
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0 

sin  n 9 

_wi_ 

0 

1 

/ 

where  /?.  is  the  rotation  of  the  ‘nodal  normal’  in  the  circumferential  direction.  The  strains 
are  now  defined  in  terms  of  global  coordinates  and  displacements  as 
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(22) 


By  a  transformation  similar  to  Equation  12,  these  strains  can  be  determined  explicitly  in 
terms  of  £,  17  and  9  coordinates,  using  the  Jacobian  matrix  [j]. 


The  strains  in  the  directions  of  the  local  set  of  axes,  one  of  which  is  normal  to  line 
V  =  constant  and  the  other  tangential  as  shown  in  Figure  7  are  obtained  by  a  transformation 
similar  to  Equation  14, 


€z  T  Yz'r' 

1  n 

t  rz'e 

€z 

1 

T^zr 

1  i 
Trz9 

2Xz'r'  €  / 

_t  K-e  2  Vs 

J -  y. 

2  Yr  9 

2  Xzr 

€ 

r 

iv 

2  Yt9 

€o 

T  ?z9 

\yr9 

€9 

(23) 
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where  again  the  matrix  consists  of  appropriate  direction  cosines  between  the  axes 

(3  x  3).  Thus  the  of  Equation  15  occurs,  together  with  unity,  because  the  circumferential 
axis  is  unchanged. 


The  vector  j  is  found  directly  from  Equation  23,  and  the  calculation  proceeds  as  before. 

Again  neglecting  the  ‘thickness*  direction,  the  [d'J  matrix  for  an  isotropic  material  now 
becomes  _  _ 
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i/ 
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0 
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0  0 
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0 

0 

0 

0 

\-v 

2 


(24) 


Clearly  the  correction  factor  k  =  1.2  must  not  be  included  in  the  last  term,  which  refers 
to  a  shear  strain  component  z'Q  in  the  plane  of  the  shell. 


SECTION  VI 

A  GENERAL  CURVED  SHELL  ELEMENT 


The  preceding  sections  describe  in  some  detail  the  process  of  degeneration,  as  it 
applies  to  basically  two  dimensional  elements.  As  mentioned  in  the  introduction,  essentially 
the  same  procedure  can  be  adopted  in  fully  three-dimensional  elements,  giving  a  general 
curved  thick  shell  element  of  arbitrary  shape  as  shown  in  Figure  8. 
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A 


Figure  8.  A  General  Curved  Shell  Element 


Once  again  the  coordinates  of  the  top  and  bottom  points  associated  with  a  node  defines 
the  coordinate  within  the  thickness  of  the  shell  and  the  directions  of  the  ‘  normal’  before 
distortion. 
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The  element  curvilinear  coordinates  are  related  to  the  global  z,  y,  z  system  by  a  relation¬ 
ship  similar  to  that  given  by  Equation  1 


X 

y 

■  I  N,  <l+2^ 

’  *i" 

yi 

f[N,  ('  "f  } 

|  2 

xi 

y  j 

z 

z . 

1 

top 

Z  j 

(25) 


bottom 


in  which  however  the  N.  are  shape  functions,  now  given  in  terms  of  £  and  17  i,e, 

Nj  =  N.  (  £  ,  tj  )  (26) 

Once  again  linear,  parabolic,  cubic  etc.  shapes  can  be  adopted  for  a  middle  surface 
and  suitable  shape  functions  generated.  These  are  as  used  in  Reference  3  for  two-dimensional 
analysis  and  are  not  given  here  explicitly. 


The  displacement  pattern  has  now  to  be  assumed,  and  it  is  slightly  more  complex  than 
one  might  expect  from  Equation  8,  as  two  rotations  of  the  *  normal’  about  axes  orthogonal 
to  it  have  to  be  determined. 


/\ 


A 


If  v„.  represents  a  unit  vector  in  the  direction  of  the  nodal  ‘normal’,  and  if  v„ .  and  v„. 
3i  li  2i 

are  unit  orthogonal  vectors,  then  we  can  write  for  the  displacements  in  global  axes  x,  y,  z 


u 

V 

11 

M 

z 

*  ui  ‘ 

V  . 

J 

w 

.  wi . 

+  I  Ni  £  ~2"  [  VH  v2il 


P: 


(27) 


where  u.,  v.  w^  are  the  displacements  at  the  midsurface  node,  and  a.,  P,  are  the  two  rotations 
about  directions  "vg.  and  -vl. .  respectively. 


An  infinity  of  systems  of  cartesian  axes  exist  which  include  the  specified  vQ.,  and  a 

Ol 

consistent  definition  giving  unique  axes  must  be  adopted.  In  general,  of  course,  the  vectors 
v^  and  Vg.  are  only  approximately  tangent  to  the  mid-surface. 


Having  defined  the  displacements  we  can  obtain  the  global  strains,  and  transform  them 
into  a  set  of  orthogonal  directions  two  of  which  lie  tangent  to  the  surface  C  =  constant.  This 
new  set  of  directions  must  not  be  confused  with  the  previous  approximate  ‘  normal’  set. 
Transformations  follow  the  pattern  of  Section  III,  but  are  not  given  here  in  detail;  the  reader 
should  consult  a  separate  paper  (Reference  11).  In  this  connection  it  should  be  remarked 
that  considerable  utilization  was  made  of  the  matrix  vector  comprehensive  subroutine  scheme. 


556 


AFFDL-TR-68-150 


SECTION  VII 

CURVED  MEMBRANE  ELEMENTS 


To  cover  the  full  range  of  shell  thickness,  and  occasionally  for  other  purposes  (when 
treating  thin  linings  attached  to  solid  elements),  one  needs  elements  which  present  no 
resistance  to  bending.  If  a  shell  is  very  thin  indeed,  the  strains  must  be  sensibly  constant 
through  its  thickness,  and  equal  to  those  at  its  midsurface. 


Immediately  it  follows  that  displacements  for  the  general  membrane  element  become 
simply 


(28) 


u 

'  uj 

V 

=  I  Nj 

v  i 

w 

.  wi 

without  the  rotations  of  the  normals  being  considered,  so  that  the  degrees  of  freedom  reduce 
to  only  three  direct  displacements  per  node;  thus  the  membrane  element  is  perfectly  isopara¬ 
metric  and  retains  all  the  properties  (Reference  3).  The  remaining  details  of  the  transforma¬ 
tion  follow  the  standard  pattern  described  in  Section  V.  Only  three  in-plane  strains,  two 
tensions  and  the  in-plane  shear,  contribute  now  to  the  stiffness  matrix. 

In  the  same  way  axisymmetric  membrane  elements  are  easily  formulated,  for  example 
by  removing  the  rotation  term  from  Equation  8. 
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SECTION  vm 

SOME  AXISYMMETRIC  EXAMPLES 

The  first  example  is  designed  to  test  the  accuracy  of  the  new  elements  when  dealing 
with  very  thick  shells.  A  simple  Lame  cylinder  under  internal  pressure  is  considered  and 
some  results  for  hoop  stresses  shown  in  Figure  9.  a.  Only  one  element  is  used  in  solution 
(and  obviously  it  does  not  matter  whether  it  is  linear,  parabolic  or  cubic).  With  extreme 
thickness  to  internal  radius  ratio  of  one  the  linear  distribution  of  stresses  fits  almost  the 
best  approximation  to  the  Lame  curve.  For  larger  radius-thickness  ratios  the  finite  element 
results  are  indistinguishable  from  the  exact  solution. 

In  Figure  9.b  a  similar  problem  of  a  thick  sphere  is  dealt  with.  Here  one  cubic  element 
is  again  used  and  it  is  remarkable  how  well  this  representation  approximates  to  the  geometry 
and  to  the  final  results. 

The  second  example  deals  with  the  classical  Timoshenko  (Reference  12)  solution  of  a 
spherical  cap.  Figure  10.  Geometrical  details  are  shown  there  and  with  a  finite  thickness 
the  distinction  made  between  loading  being  applied  to  outer  (Case  I)  or  inner  (Case  II) 
surfaces.  The  average  stress  resultants  reproduce  almost  exactly  the  theoretical  curve  using 
a  division  of  24  cubic-type  elements,  but  it  is  of  interest  to  note  the  differences  arising 
between  the  two  loading  cases. 

Figure  11  shows  the  same  problem  solved  for  a  thicker  shell  (t  =  9").  Results  are 
compared  with  a  full  two-dimensional  solution  of  type  indicated  in  References  3  and  13.  The 
computer  times  necessary  for  equivalent  solutions  here  were  in  the  ratio  of  1  :  3. 

The  third  example  deals  with  a  very  thin  restrained  pipe  (Figure  12)  which  has  also  been 
previously  solved  by  finite  element  methods  by  Grafton  and  Strome  (Reference  8)  and  Percy 
et  al  (Reference  9).  Here  the  effect  of  various  subdivisions  into  cubic  elements  is  studied. 
Even  with  two  elements  the  displacement  pattern  is  well  reproduced  but  six  elements  were 
needed  to  obtain  reasonable  stress  predictions.  For  solution  of  the  same  problem  in  Ref¬ 
erence  9,  a  much  larger  number  of  elements  was  needed  even  to  obtain  displacement  patterns. 

The  last  example  concerns  another  test  example  now  of  an  axisymmetrie  tube  under 
nonsymmetric,  harmonic  loads.  Here  the  problem  and  the  results  for  several  harmonics  as 
given  in  Figure  13  and  are  compared  with  results  of  Percy  et  al  (Reference  9)  as  well  as 
those  of  Budiansky  and  Radkowski  (Reference  14). 
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Figure  9.  Lame'  Solution  for  Hoof  Stresses  (a)  Cylinder  Under  Uniform  Internal 
Pressure  (b)  Sphere  Under  Uniform  Internal  Pressure 
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Figure  13.  Cylinder  Under  Harmonic  Edge  Couples.  Finite  Element  and  Numerical 
Results  of  References  9  and  14  (Solid  Line)  Compared  With  Present 
Solution.  Analyzed  Half  With  15  Elements  in  A.P.  (1st  Element  is 
0.02  From  the  End).  =  Meridional  Bending  Moment  (Amplitude) 
u  =  Radial  Deflection  (Amplitude),  E  =  1,  v=  0.3. 
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The  general  agreement  of  displacements  and  moments  is  excellent  with  exception  of 
the  displacement  results  for  the  fifth  harmonic.  This  discrepancy  remains  somewhat  puzzling 
since  a  full  two-dimensional  analysis  (Reference  13)  of  the  same  problem  also  differs  at 
the  same  point. 


SECTION  IX 

SOME  ARBITRARY  THICK  SHELL  SOLUTIONS 

The  first  problem  is  that  of  a  circular  ring  line  supported  at  the  base  under  the  action 
of  gravity.  Figure  14.  Results  agree  well  with  those  predicted  by  simple  engineers’  bending 
theory. 

A  more  complex  situation  is  illustrated  for  a  variable  thickness  tube,  supported  at  two 
points  as  in  Figures  15  and  16.  Again  the  loading  is  that  due  to  self  weight. 

The  last  example  is  that  of  a  fairly  thick  shell  representing  an  arch  dam  of  symmetric 
shape.  Figure  17  shows  the  general  division  into  nine  thick  shell  elements. 

In  Figure  18  the  downstream  displacement  plot  is  compared  with  results  of  a  complete 
three-dimensional  solution  of  this  problem  (Reference  4).  The  results  agree  very  well 
though  it  must  be  mentioned  that  a  slight  difference  in  the  manner  of  load  application  existed 
between  the  two  solutions. 
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Figure  15.  A  Variable  Thickness  Conduit  Supported  at  Two  Points  Under  Action  of 
Gravity.  Analyzed  With  32  Parabolic  Elements  (4x8  Mesh). 

(a)  Geometry  (b)  Variation  of  Displacements  w/Xd4  and  Axial 
Bending  Moments  M^/Xd  Along  Top  and  Bottom  Axes. 
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Figure  18.  Downstream  Displacements  of  Air  Face  on  Crown  Section  of  Dam. 

Comparison  of  Present  Solution  With  Three  Dimensional  and  Thin 
Shell  Solutions  of  References  4  and  15. 
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It  is  of  interest  also  to  show  a  comparable  solution  based  on  triangular  finite  elements 
and  classical  thin  shell  assumptions  (Reference  15).  The  difference  of  results  due  to  effects 
of  shear  deformation  now  included  is  significant. 

The  stress  agreement  (again  for  the  crown  section)  is  indicated  in  Figure  19. 


SECTION  X 

CONCLUDING  REMARKS 


The  new  element  has  been  shown  to  be  capable  of  dealing  adequately  with  thick  as  well 
as  thin  shells.  By  virtue  of  its  curved  shape  it  can  fit  any  geometry  adequately. 

At  one  end  of  the  scale  it  duplicates  thin  shell  solutions  and  is  of  comparable  economy 
with  other  elements  used  in  this  context.  At  the  other  end  of  the  scale,  it  is  capable  of  dealing 
with  very  thick  shell  situations  much  more  economically  than  by  corresponding  two  or  three 
dimensional  formulations.  Its  particular  use  is  envisaged  in  the  complex  practical  situations 
where  thick  and  thin  shell  elements  are  combined  with  possible  membrane  regions. 

The  new  element  is  being  incorporated  in  the  NIFE  (Numerically  Integrated  Finite 
Elements)  system  developed  at  University  of  Wales,  Swansea. 
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